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INT R O DHOTI OH 

1,1 Introductiaa: 

Plow induced vibration is a wide field associated xd.th. 
engineering structure placed in flowing fluid. Energy needed 
for this type of vibration is derived from flowing fluid around 
it. The vibration may be caused by far field noise or near 
field noise. In the first case turbulence, eddies and cavitation 
etc, plaj?' majoi'* part ^iiile in the second case it is mainly 
due to pressure fluctuation in boundary layer of the submerged 
body. . ' 

Plow induced vibration problems were studied earlier for 
fatigue damage of nuclear reactor core where the flow velocity 
is comparitively high for hi^er rate of heat transfer . Plow 
may be parallel or perpendicular to the axis of the rod, Por 
cross flow which is mainly important in heat exchanger, the 
vibration is due to the far field noise. The problem has been 
widely studied by several workers and empirical relations are 
available as design guide. However, in the case of parallel 
flow B-ome theoretical models have been attempted using statis- 
tical me'bhods. Due to lack of infoimiation about the flow 
field around the rod, it is assumed that the rod is excited 
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by a random time varying pressure field. Statistical correla- 
tions have been derived to obtain a qualitative idea in this 
field. 

Ihe present work is attempted to study the vi-braticn of 
long nuclear fueld rods (PVffi reactor) which are supported 
periodically by grid, plates. Effect of continuous eirosion 
of baffles resulting in enlargement of baffle holes and notching 
of fu.,.l rod due to hammering and sawjtng of the vibrating rods 
are observed. A study of the problem requires the Icnowledge of 
i»he response of a periodic structure unde^r random convected 
pressure loading, 

^•2 Revie w of Previous ¥ork si 

Ihe interest in the area of flow induced vibration was 
reactivated by Ashley and Haviland [l] in connection with 
$rans Arabian pipe lines. A lot of theoretical work in this 
area has been taken up since then. A nice comprehensive 
review of the literature is given by Paidousis and Issid [2] . 

A brief review of the e^^perimental work is given below. 

Burgreen et al. [33 studied the vibration of rods with different 
end conditions in parallel flow and it was observed that the 
vibrabicn is of self excited t 3 rpe , The rods vibrate at their 
own natural frequencies independent of flow velocity, Quinn [4] 
also measured the vibraticai of a submerged rod in water at 
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temperature ranging from 21 to 61°C* The experiments were 
also carried out for two phase flow. Paidousis [5] performed 
some experiments on a single I'od and a bunddle of rods in 
paral3-el flow. Based on his experiment he revised the empirical 
relation previously derived by Burgreen et al. Reavis [6] in 
1969 offered an attractive approach to the problem postulating 
that the vibration is excited mainly by boundary la,yer noise. 

He compared his results with the theoretical, solution given by 
Thomson [7] for the first mode response of a pinned beam to a 
distributed random loading. The forcing function as given by 
Bakewell [8] was used in the previous work, Balcewell suggested 
that the pressure fluctuation is convected down stream and is 
a function of frequency. Ram Jiyavan [9] tried to derived 
some statistical correlations of the pressure characteristics 
of parallel flow over a, thin rod. He also observed the spectral 
variaticn of pressure field \ri.th flow velocities, ill these 
works that have been done in the field of parallel flow induced 
vibration are on single span i*od or tube. But when the nuclear 
fuel element is supported intermediately the problem is modelled 
as a periodic structure under random distributed loading* 

The most common method fcr analysing a structure subjected 
to random loading, is the modal method. In this method, it is 
assumed that the modes are independent of each other -whibhci 
is not true for a periodic structure. Moreover, so many modes 
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are associated -with periodic structure with, many supports 
that the method becomes unx^ieldy. Lin [10] applied tx'ansfer 
matrix method to solve the problem of beam xm.th many supports. 
Ihe use of finite element method far solving such problems 
has also been attempted [11] . Recent application of wave 
propagation theory in periodic structure by Mead [12] made 
the problem easier for analysis. The main advantage of this 
method is its simple analysis and the negligible computational 
time wlxLch is independent of number of supports. S.en Gupta [15] 
developed a graphical Viethod, using wave approach to determine 
the natural frequencies of a periodic beam. Mead [14] used 
the wave approach to determine the response of a periodic 
structure xmder convected harmonic loading. He shox^ed that 
the ’coincidence’ occxnrs at a frequency when the convection 
velocity of the pressure is equal to the phase velocity of the 
free waves propagating -the. structure . The limitation of the 
theory is that a closed form solution is possible only when 
the beam is uniform. Hoxrever, Mead and Mallik [15] have 
developed an approximate method using ’assumed mode’ tedmique 
to determine the response of an infinite periodic structure 
with non-xmiform element. Rao [16] applied this method to 
find the response of a finite periodic beam. 

1.5 Scone of Present Work ; 


The present work deals with the vibration response of a 
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long i-od which, is supported at regular intervals. Auto 
correlation functions and power spectra of both pressure 
(excitation) and vibration (response) have been obtained. 

.¥ith these results frequency response functions of the system 
have been computed with different flow velocities. Results 
have been obtained in various spans of the rod. The experi- 
mentally obtained results are compared wi,th the theoretical 
pi*edictions obtained by using imve approach. Auto correlation 
coefficient of both pressure and response have been computed 
using Honeywell SAT— 48 correlator. Poi-jer spectral density 
functions have been comiDuted from auto correlation results on 
DBC-10 using the well known relation. 

P(f) = Lim i c(t) dr (1.1) J 

X ■><» 


where c(t) is the auto correlation coefficient and f denotes 
the frequency in Hz, 

¥ith H discrete values of c(t) spaced at A t apartj 
the poTOr spectral density function can be expressed as 


P(f) 




c(o) + 


M 

2 Y 

^ c(o A t) cos 2%f2 A t * 
j=l 

M 


At c(o) + 2 At ^ c(j/^t) cos 27cfjAt 

( 1 . 2 ) 
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Haring power spectra of both eKoitatlon (pressure) and 

response (vibration), the system frequency response function 
can be computed. tie relation 


S^,(f) = |H(f)|^ Sj.j,(f) 


(1-5) 


where 


spectrum of response, 
power spectrum of excitation. 



mimimmjih set up 


2.1 Introduction: 

A sdaematic view of tlie esper imental set up is shown in 
Fig, (2,1), The horizontal test loop consists of several 
remova,hle sections of 15 •5 oni dia C.I, pipe, A centrifugal 
pump (22 H,P.) delivers wa,ter to the test loop. Flow of the 
water can. he controlled hy three valves (¥^, V 2 * * 

outlet water is recirculated to the reservoir again. 

The test loop has been designed for minimum far field 
noise. The loop consists of two flow straightners , one AFT 
filter and two sets of acoustical filters (F^j T’2^ • 
acoustical filters are placed in upstream and down stream of 
the test section to take care of up stream noise and reflected 
noise as well. The outlet water passes over a V notch while 
flowing through the re circulating- channel , Detailed discu- 
ssion of the test loop is given in ref. [9]. 

2.2 Test Section; 

The test section consists of a 2m long, 15.5 cm, I.D. 

M.S, pipe placed between filters F^ and F 2 * It contains a 
^5 ft, long test elan ent periodically supported from the outer 
pipe by pointed supports. The supports are designed to simulate 
the simple support condition as nearly as possibie « The test 
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element is push fitted with wedge shaped rings at 1 ft. interval 
to provide a line suppox’t on the test element, lach ring is 
supported by four pointed end bolts so that the test element 
rests coaxially t-ri.th outer pipe, Ihe test element is made of 
five parts of equal length. Bach part is connected to its 
adjacent part by means of threaded socket, The support locations 
are so chosen that the sockets are at middle of each bajr. 

The sockets are instrumented with strain gauges to measure 
the vibration response at the centi*e of each bay. The above 
arrangement has been done for some specific advantages. 

(1) It is easy to instrument the small socket rather than the 
long test element, 

(2) One instrumented socket can be changed from one bay to 
another so that we need not instrument all the bays at 
a time , 

2.3 Instrumentation and Signal Conditioning: 

The experimental procedure consists of two parts: 

(a) Measurement of pressure which is the excitation 
function. 

(b) Measurement of vibration response , 

Bor measurement of pressure, two pressure transducers (BMDBYOO 
model 2520) have been moxmted on the walls of the outer; pipe 
with the assumption that the same pressure fluctuation is felt 
on the pipe wall as on the test element. Signals from pressure 



transducers are fed to the charge amplifiers, Prom charge 
amplifier the sigaal. is taken to the correlator via a low 
pass filter. 

Vibration measurement has been done with strain gauges. 
Since the amplitude of vibration is very low semi-conductor 
strain gauges are used due to their laigh sensitivity and gauge 
factor. Two strain gauges have been mounted on the top and 
bottom of each socket, leads of the strain gauges have been 
taken out throu^i the hollow test element so that the flow 
over the test element is not disturbed, leads of the strain 
gauges are taken to a Budd strain indicator (Model P~350) to 
balance a half bridge circuit. The output from the strain 
indicator is fed to a linear amplifier in order to amplify 
the signal. This amplified signal is fed to the correlator 
via a low pass filter, 

2,4 Bsperimenta l Procedure ; 

2,41 Sampling; 

Since SiI-48 correlator is a digital unit, choosing 
proper sampling time is important. Two t37pes of difficulties 
may arise for improper choosing 'Of sampling tiii©. If the 
sampling time is too short, to obtain the necessary information 
from the data will involve too much time and increase the cost 
of calculation. On the other hand sampling at points which 
are too far apart will lead to statistically poor data :: 
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especially for low and high freciuency components of the original 
dara. This is called aliasing. In general sampling time is 
chosen from the knowledge of highest freq,uency present in origi- 
nal record. The relation hetween highest frequency and sampling 
time is given hy f^ = where f^ is called lyquist folding 
frequency and h is sample time interval. Prequency above 
hyquist frequency T'ri.U. be folded back. In our experiment yiH 
che data have been taken below 100 cps. Therefore the proper 
sampling time to avoid aliasing is 5 ms. To be on the safe 
side so that no frequency component is lost the sampling time 
in most of the re cording is taken to be 1 ms . Si0;-48 Correlator 

has a series of sampling time ranging from '1 ps to 2 aes and 
gives an output of 400 correlated data points. 

2.42 Measurement of Correlation functions: 

As mentioned earlier that the experimental procedure 
consists of two parts: (1) Measurement of prepare (2) Measure- 
meno of vibration resj)onse. Por eadh part auto correlation 
functions have been computed in SAC -48 correlator. aAE-48 
performs correlation for 400 points , Signals from transducers 
are fed to either channel of the correlator. Appropriate delay 
is introduced in one of the channels for computing the correla- 
tion function. Digital output from correlator are taken to a 
digi t/Sl computer (DEC~10) for computing power spectral density 
funcbion, Por measurement of convection velocity cross corre- 
lation functions have been computed between two pressure trans- 
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ducsrs placed axially apart. Cross correlation function 
aceurately deterraines the delay time {T^) of the signal 
"between txTO pressure transducers. Convection velocity cai 
he found from relation Cp = ^ , \&ie:ce Gp denotes convection 
velocity and d represents transducers separation distance. 

The cut off frequency in most of the measurement is kept helow 
100 Hz because the strain gauges are unable to respond at 
higher frequencies. All these experiments have been done for 
different flow velocities. The flow velocity has been measured 
conventionalljr from the head of the water over a T notch while 
it was passing through the recirculating channel. 



TIEBORETICil. MiilSIS 

3.1 Introduction; 

A poriodic s cruoture consists of identical elements connected 
in iaenuical manner* A t 3 ^Dical example of this structure is a 
multi supported continuous heaii: mth identical supports at 
regular interval, She response of such structure was analysed 
in the pasc "by modal memhod* Vi/heii the beam contains many supports^ 
so many modes are to he included_^he modal analysis becomes 
cumbersome. Under such circumstances it is convenient to consider 
the structure as infinite so that discrete modes can be ignored. 
Ihe motion can be conveniently analysed by a so called wave 
approach. 

3.2 Pundamentals of ¥ave Approach; 

A simple wave propagating along an one dimensional damped 
s cructure can be characterised by the wave number K = 2%/\t 
wlaich is nothing but the phase difference between wave motion h 
at two points which are unit distance apart, fhe wave decays 
as it propagates. The amplitude ratio of two points which are 
unit distance apart is e*^. The quantity (6 + iE) is called 
'che propagation constant which defines the amplitude and the 
phase relation per unit length of the wave motion. Suppose 
an ixifini ce beam with equispaced supports is excited at a 
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single point P "by a liarrnorD. c force The wave motion 

>ri.ll propagate outwards fx’om P. The displacements in any two 
bays imjT- not be identical. Hoviever, the amplitude ratio and 
pha.se difference of two identical points in ad 3 a, cent bays are 
same, represents the phase difference between two identical 
points in ad3‘acent bays over a distance S. while e represents 
the amplitude I’atio . 

The characteristic propagatiai constant is now +(112, + iPj^)« 
The negative sign implies that the wave is decaying as it is 
propagating from left to right while the positive sign implies 
the reverse direction. 


3.3 Properties of Propagation Constant: 


Consider two adjacent bays A-B and C~D resting on rigid 
simple supports, restrained at each end b5?' a rotational spring 
(figure 3 . 1 b). The total rotatimal stiffness at each 
support is Suppose the beam is subjected to free vibration. 

The moments at (r~l),r, (r+l)th supports are 
respectively. These are all harmonically varying complex 
quantities, let us consider the span A~B, slope of the beam 
at point B in span A-B may be expressed as 



M 

B A ^ 




(3.1a) 


where p is called receptance function which is defined as 
follows 




FIG.31 PERIODIC BEAM ELEMENT 
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^BA ~ point B due to unj.t nionent at point A. 

^BB ~ slope at 3 due to unj.t moment at point B. 

Similarljr, if ire consider the span C-D, the slope at point C 


~ ^CC \ ^OD ^'Wl C3.1h) 

If the beam is uniform and resting on identical supports 
then, Pt 33 = -Pqq and = "*^CD* 


Bor caitinuity of slope = 0^. Hence, equating equations 
(3,1a) and (3.1b) and using above receptance relations, ¥e get 


M 4. f .SS Jvf .1, J/j 

Pba ^ 


0 


(3.2) 


The above relation applies between two adjacent free bays. 
Such set of equations satisfy recurrence relations like 


H. 


r+l 




IL 



(3.3)^ 


Substituting these rela,tions in equation (3.2), ^f7e g-et 


u 

e‘ + e 


1-1 


■ 2p 


or 


Goshhu 


SB 
^BA 


Ebb 

’Pba 


( 3 . 4 ) 


Ihe right hand side of the equation is highly frequency 
dependant and complex in nature if damping is present in tbe 
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P "D p 

beam. Por an undanped beam^ if ~ --■y.-i is greater than 4-1, p. 

^BA 

is 9ntirel3r j^eal. Therefore, the wavs is decaying in nature 
but there is no phase difference between two adjacent bays* 
The wave is non propagating t57pe. If ~ is less than ~1, 

pba 

[.I is of the form p. = ja.^^-rin:. The waves in two adjacent bays 

p '~>" Q 

are in counter phase but still deca 5 nng in nature. If - 

*^]3A 

is between -1 to +1 then p, is entirely imaginary. This means 
there is no decay in wavs amplitude and wave can propagate* 

P -p-p 

The frequency range in wliich - lies between -1 and +1 is 

Pba 

called the propagation band. 


3.4 Res^Donse of An Infinite Periodic Beam under Harmonic 
Convected loading; 


Consider an Euler -Bernoulli beam resting on periodic 
suppoi'ts ai distance Si apart and excited by a harmonic pressure 
field such that the force per unit length of the bean is 

P(x,t) = 


where K is the wave number and is equal to ^ with Cp as the 
convection velocity. Phase difference between two points at a 
distance x apart is -Kx. 

The motion of the beam is governed by the following equation 


BI — -j* + m ^ 
gx^ 3 1'^ 


^l(aj.t-Ex) 


(3.5) 


If damping is present in the beam, then El should be replaced 
by BI(l+ir]) ' viiere p = loss factor. 
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BI _ flexural rigidity of the beam eind m = mass per unit 
leng-hi of the beam. 

let us put =c = C<l/2, where ? Is the non-dlmensloi^l x co- 
ordinate, with the origin at the middle of the bay. 

Assuming the steady state solution W = W (j)e^'*' and 
subtltuting this in the governing equation one gets 


16(SlA‘^)(l+ir,) af| 


dC 


2 ~ ■ ~iE£C/2 

m <»> y = P e 
0 


or 


d^y ' ^ S,^ W 


^ a 
o ^ 


-iEilC/2 


dC 16 BI(l+ir)) 16 Bld+ip) 


(3.6) 


(3.7) 


let ue introduce the foil o.ring non-dimensional quantities 


( a) frequency parameter fj - 

SI 

(b) loading parameter P 


!o5 

‘bi’ 


( c) convection velocity . parameter CV = Cpji 
Substituting these quantities in equation (3.7), we get 

n 


dZl 


n2 - 

n ifif 


P Ji 


£^^4 16(l+iTi) I6(l+ir)) 

The solution of the equation is 


-ipQY 5 


(3.8) 


n=l 


(3.9) 
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^ ( l+XT) ) 2 


where Pjj, is taken -bo be 1. 


are the four fourth roots of the quantity ^ . 

A^are the constants to be found out from the boundLy conditions 
llov consider the notion of a :free bey. Since there is no loading 

in free hay the right hand side of the equation 3.8 is absent. 

Hie solu'fcioii ill sucli case is 


n=l 


>n € 


(3.10) 


To find the receptance functions p's, apply a unit non dimen- 
sional harmonic moment at end A of the span IB figure (3.1c). 

Ihe boundary conditions are: 


t (-1) 
I? (+1) 


■W’'(-l) 


~i) 


( l+ip ) ( ■2.,. { l+iT] ' 


(3.11) 


■¥’ »(+l): 


IK t’(+i) 


2 ( l+ip ) 


Here we have considered that the supports are transversely rigid. 
Ine quantity is non dimensional rotational stiffness at eac±L 


support 
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Q-oov'e ooundary conditioiLS can. lae solved uiiiq_uel 3 ?'» 
HoM the re cep -bailee and are given by 


4 

I 

n=3- 


?AA = »’(-!) = I \ An " " = -S 


BB 


Pba = ^ \ h ^ “ = -P 


( 3 . 12 ) 


n=l 


'iB 


In general, ^j^’s are complex and frequency depoident, (Therefore, 
receptance functions are also complex and frequency dependent. 
Using these values of receptance function propagation constant 
(p) can be calculated from equation (3.4). 

3.5 Borced Wave Response in An Inflni-te Beam: 


In order to find the forced i-rave response generated in an 
infinite beam, let us consid.er the equation (3*9) where P.I. 
is not zero together with 1iie boundary conditions given below. 

¥ (- 1 ) = 0 

W (+1) =0 (3.15) 

Slopes a"t either end should satisfy ¥’(+1) = e ¥’(-1) 
and bending moment rela'bion is given by 

El (1+ip) ¥ (+1) + I ¥ (+1) = e El(l+ip)¥' »(-l) 

Using above boundary condition A^s can be determined and the 
displacement is given by equation (3*9). 
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3 


Response of A Unit 
Held: 


g Periodic Beam, to Coiivected Pressare 


Pile response of a finite periodic beam consists of tTO 
parts: 


(a) Forced wave notion generated in an infinite periodic beam. 

(b) Response oori'espondlng to reflection of waves from extreme 

6iids wliicii arc0 ©33 8111/1 Elly fr©© waves# 

Til© solution of part (a) is given in the previous section# 
xhe solution ox part (b) for free waves can be written as 


4 ^ 

Wf (5 ) = I 
n=l 


B^’s are the constants to be solved from the foHox^ing 
boundarj?- conditions. 


(- 1 ) =0 
(+ 1 ) =0 

Slopes at either end is related by the relation 

¥»C+1} = ¥»{-l) 


The above three boundary conditions can not determine four 

Bo 

bT = 1*-X, 
. Bt 


B. 


constants uniquely. However, if we consider ^ 

B« 1 

% 


bo f 


then B^*s are determined uniquely 


0 "1 

apart from a common factor i^iicli comes in the am- 

* 





wave 



The above bouMary conditione caia be re-written ^ the form 

4 ~X -X 

T b e ^ = -e 

^ n ■ 

n=2 

4 ),n *1 

‘ h® = -9 (3,14) 

n=2 

IJ’ 1 

Z \ - ® ® > = - V® - " •® ^ 

n=2 

Prom these equations -with given and -'for positive and 
negative going waves t-wo sets of values b^_^ and can be 
;found • 

Total response of a beam is gi-ven by 

WjCe) = V (5) + yf+(5) + WfJ?) 

Where wCj) is the forced response in an Infinit® 163111 

X 5 X 

¥^(e ) = ¥ ) + B^ l e ^ + ^2 I ® ^ ^ 

(3.15) 

The unknown constant and B2 can be calculated from exureme 
boundary conditions. For a five span beam, these are as follows 

Sl(l+ir|) ¥» (-1) = (-1) 

-BI(l+iTi) ¥” (+1)^“ ¥J {+1)5 
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Substituting the value of Wj from equation (3.15). w get 
j3I(l+ir|) -l)+WJ^(-W+''fJ-l>] = 

= E^[e'5^^"w-C-1) + eS'^t W^^(-l) + e^l^- ¥^_(-l)]. 

i.-^o -n"^ nni^ tC can be determined. It should 
i^om above equaticns and nj_ can ue u 

be remembered that the response given by eq,uation (3* 

for the 1st bay only. 

The response in the nth bay is given by 

^ ^-i(n-l)Kl ^( 5 ) ^ Wf^(C) t 

(3. IT) 


f 1 for iiio *0111 "t valuo 

The response given by equation (3.W) is lor 

Of the loading parameter i.e. P = 1- *‘>1^ 

consider 1“ non-dimensional form (i.e. toviduig 

by the span length l). Ibe second derivative of ,.( 5 ) 
the non-dimensional curvature (bending strain or stress . The 
amplitude of this quantity is nothing but the frequency response 
function 1 H(f) 1 if we consider the bending strain as the ou pu . 
In the present work, this lH(f)| is calculated for yrio-^ 
values of OT at the .mid-point of the spans (I.e. - 

both theoretically and experimentally. 



QHAPTBR - 1.1 


EJ3SULTS MD DISCUSSIONS 

The test element used in the present work has the 
following physical and geometrical characteristics: 

Outer diameter = 23 mia. Material - Mild Ste 

Inner diameter = 20 mm. 1st natural frequency - S1& m 

NO. Of span =5 Modulus of elasticity, E = 210 

Span laigth = 304*8 mm. OM/m 

agures (4.1, 4.2) the auto conrelatlou functions of 

pressure and vibration for different flow velocities. The 
curves decay with increasing delay tine showing the randomess 

in the process. The curve become steep with mcreasuig velooi y. 

4 . a+T on curve is an indication 

Since the slope of the auto correlation curv 

. . o-nni it shows that highsr 

of frequency components present in a g » 

frequency components are introduced at higher velo 

amplitudes of the curves change their si^ but the 

the curve is positive satisfying the condition of positive 

definiteness. 

Pigure (4.3) bhows the PSD curves of procure signal. 

oitinntsd in a digital computer using auto 

These results are computed in 

correlation data with the help of equation (1.2). 
that most Of the energy is associated wi«. 1» 

components. The curves *ov peaks between 30 and 
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Li.'iinly due go ecldiss ond. 'fcux’lDuJ.eiice « Tiie psa,lcs a.re 
siiiiT’feiiio to ~Qhe liighei* frequencies as tlie flow velocity 
ill ere uses « Hii. s sIiowg tlie flow dependence of th.e pressure 
field. At higher flow velocity the eddy life time decreases 
there by showin-o; shift in the frequency spectrum. 

Figure (4 •4) shows PSD curves of the vibration response. 
Phese results are also computed in the same manner as the 
pressure spectrma, Phese curves also show peaks between 30 
and 40 Hs. Phis shows the turbulence and eddies are primarily 
responsible for enciting the rod. 

Figure 4»5 ^ows cross correlation curves between two 
pressure transducers which are placed axially apart. Phe axial 
separation between the transducers is denoted by (d) wMle 
P^ represents the delay time corresponding to peak amplitude 
of the cross correlation curve. Phe convection velocity 
is calculated from the simple relation, Cp = It is found 

that when the flow velocity is 49*6 m/sec the convection 
velocity is 46 aa/sec. Phis shows that the disturbance in the 
pressure field is travelling with a velocity almost equal to 
flow velocity. 

Figures (4.6, 4.7, 4.8, 4.9) show the frequency response 
function H(f2) of the system. Phese results are computed from 
PSD data of pressure excitation and vibration response using 
equation (1,3). Phe carves are plotted against non~dimensional 
frequency a(= aicfJl ^ for different non-dimensional 
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convoctioii volocitieg 07 (= C# VkT \ a- . . 

r gj / • Since the absolute 

measuuenient of the amplitude was not ‘done, the experimental 


curves are compared mth the theoretical curves after matching 
the peaJ: amplitude value. It is seen that iiie relative 
magnitudes of pealc and valley are almost same. On ccmparing 
the experimental results with those obtained from theory, it is 
found that there is a shift in peak frequency between the 
experiraontal curve and the theoretical one. Ihe experimental 
curves show pealcs at higher frequencies than the theoretical 
ones. Ihis discrepancy may be due to the fact that the 
pressure has been measured on the wall of the outer pipe 
rather than on the surface of the test element. In such cases 
the measured pressure may not carry the true spectral informa- 
tion of the actual pressure* Moreover pressure field disturbance 
dCEC'tocoDupling between the acting prepare and the resulting 
vibration of the element has not been accounted in the theory* 
Ihis effect seems to be quite important and produces a signi- 
ficant diange in the pressure field across the floxf. 

However to clear this doubt a pressure field record on 
the surface of a 1" dla rod was taken and it has been found 
that the frequency response curve, using the above result 
matches better with the theoretical prediction (ELgure 4*10), 

It should also be noted that in theory the supports of the 
test element are considered as simple line supports where as 
in practice some rotational stiffr.ess is always present at the 
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supports, fills d.isci'‘3paiicy also taiids to push the peak 
frequency to the higher value, 

fhe p sales of the eiiperi mental curves shift to higher 
frequency mtli increasing convection velocity which is also 
predicted 'by the theory, fhe value of the loss loss factor (ri), 
used to obtained the theoretical curves, has been computed by a 
trials and error approach, first an experimentally obtained 
curve of the frequency response function is matched with a 
theoretical curve for a particular value of damping, faking 
this value into account other results have been predicted. 

fhe frequency response curves show that there is not 
much difference between the response of the third span and 
the fifth span in 1he convection velocity range 0,05 to 0 . 05 * 
fhis is also predicted by theory. But it is expected that 
there will be a mark difference between the response of the 
third bay and fifth if the ccsivection velocity increases 
beyond a n on-dime nsicaaal value 2, 
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FIG- r.-6 FREQUENCY RESPONSE FUNCTION AT THE 
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FIG- 4-10 frequency RESPONSE FUNCTION AT THE 
CENTRE OF THIRD BAY AFTER USING THE 
PRESSURE FIELD ON SURFACE OF TEST 
ELEMENT - - 



CHAPTER - V 


GCWCEQSIOHS 

In i±Le light of the above results the following conclu- 
sions can be dr ami. 

(i) Since the PSD and the frequency response curves are 
changing with the flow conditions , these can be considered 
as system representation rather than extraneous noise. 

(ii) In the range of parameters of the present work, tlae 
vibration of the rod is forced type and highly flow 
dependent, 

(iii) A periodically supported rod subjected to parallel flow 

high 

may vibrate withj^amplitude at a frequency much lower 
than its natural frequency. This peak frequency increases 
with increasing flow velocity. 

(iv) The pressure field is mainly dependent on -eddies which 
propagate with the flow velocity. 

(v) Pressure field should be measured on the specimen even 
at this low flow velocity. 



iPPSNDIX 


Honeywell Correlator and Signal inalyser (Model Sil~48) 
Product Description ; 

Sj0;~48 is a high speed processing instrument for on 
line correlation, signal averaging and prohabiUty computation. 
Phe outstanding features of this instruments are; 

(1) 400 point analysis 

(2) 5 MHz (0.2 psec) sampling rate 

(5) Digital exponential averaging 

(4) 800 points precomp delay 

(5) Distal solecticn and read out of any bin. 

Operatio n ; 

( a) Correia ti on : 

Correlation is a fundamental tool for time domain analysis 
341-48 provides auto and cross correlation functicn with incre- 
mental time lag ranging from 0.2 psec to 2 sec. Por both auto 
and cross the processing of signal is identical. The only 
difference occurs in input circuit. In cross correlation two 
different signals are applied to the two different channels of 
341-48. In auto correlation one input signal is applied to 
either channel and internally routed to the 2nd channel to be 
automatically processed as in Cross Correlation. Eegarding of 
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whether the mode is auto or cross tho operation is as follows: 

input signal is sampled at the rate indicated hy 
chosen sample increment switch and Q.uantized in Tooth the channels* 
Channel I has got 400 word memory I'diile channel II has got one 
word memory (the present sample)* At time called present time, 
che cooorelacor performs the multiplication of present word of 
channel II with each 400 word of channel I. Ihe 400 results 
loecome the first entry of R( 0 ), H(l) . . . . . .R( 399 ) * R( 399 ) 

stands for the correlation value at log time 399 t. The operation 
is repeated for repeatedly updated values in channel II and 
corresponding updated values in channel I and then averaged. 

For other details of operations like 'Enhancement, Pro balDility 
density, Probability distribution' see Ref. [I 7 ] . 

Application : 

In this experiment correlator has bean used only for 
specific purpose like measurement of auto correlation and 
cross correlation of pressure and vibration sigial of a rod 
in parallel flow. The cross correlation results gives the 
information about velocity while the auto correlation results 
are used to compute the PSD function. Honeywell has got a 
product *SAI-470 Eourier Transfoimier ' win. di can be connected 
to SaI ~48 Gorrelaior by interface connection to get the PSD 
function directly by Fourier analysis of correlation values. 
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